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Abstract 

Ferromagnetic compound PrRu2Si2 exhibits a giant magnetocrystalline 
anisotropy of about 400 T. Its ordered moment below T C =14K reaches 2.7 
\ib and is parallel to [0 1] crystalline direction. We have attributed the 
magnetism of PrRu2Si2 to the Pr ions and performed calculations of the 
fine electronic structure of the Pr ion in the tetragonal symmetry, relevant 
to PrRu2Si2 taking into account crystal-field and inter-site, spin-dependent 
exchange interactions. Spin-dependent interactions have been taken into 
account by means of molecular-field approximation. The derived energy 
level scheme is associated with the removal of the degeneracy of the lowest 
multiplet given by Hund's rules, 3 R"4. Magnetic and electronic properties 
resulting from this fine structure are compared with all known experimental 
results. Our calculations reproduce well the zero-temperature moment, 
temperature dependence of the magnetic susceptibility, single-crystalline 
magnetization curves with the anisotropy field of 400T, the specific heat 
with the sharp peak at T c as well as inelastic-neutron-scattering data. 



1 Introduction 

Ternary compounds RM2X2 (R= rare earth M= 3d, 4d or 5d transition metal 
and X=Si or Ge) exhibit large anisotropic magnetic properties. PrRu2Si2 ex- 
hibits the largest measured magnetocrystalline anisotropy. The anisotropy field 
at 4.2 K is estimated to be so enormous as 400 T [p]]. PrRu2Si2 crystallizes in the 
tetragonal structure ThCr 2 Si2-type and orders ferromagnetically below T<7=14K 
[0. Its magnetic susceptibility shows strongly anisotropic behavior in tempera- 
tures above Tc ||. The temperature dependence of the specific heat of PrRu2Si2 
shows the well-defined anomaly at the ferromagnetic ordering temperature. The 
temperature dependence of the magnetic entropy computed from the magnetic 
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contribution shows value Rln2 at about 30K, indicating the existence of two 
closely-lying levels. Inelastic-neutron-scattering experiments [| reveal local 
excitations with energies of 2.25 and 29 meV providing strong argument for the 
existence of the localized CEF-like levels. 

The aim of this paper is the evaluation, on basis of all known experimental 
results, of the energy level scheme of the fine electronic structure of the Pr 3+ ion 
in PrRu2Si2 associated with the configuration 4/ 2 . This structure provides con- 
sistent description of main magnetic and electronic properties of the compound 
under study, the direction and the ordered moment value as well as the sharp 
peak at Tc in the temperature dependence of the heat capacity, in particular. 

2 Outline of theory 

In the individualized-electron model / electrons keep their individuality also being 
placed into intermetallic compounds. In the intermetallic compound there coex- 
ists a few (here only 2 are specified) physically important electronic subsystems: / 
electronic subsystem(s) and conduction-electron (c-e) subsystem. These two sub- 
systems are described by completely different theoretical approaches: localized- 
electron and itinerant-electron models. The / electron subsystem exhibits the 
discrete energy spectrum associated with bound states of the electronic system 
f n . Itinerant electrons occupy the conduction band states. We have attributed 
the magnetic properties of PrRu 2 Si 2 to be predominantly due to the 4/ 2 electronic 
system of the Pr 3+ ions because the c-e susceptibility is small and largely tem- 
perature independent as we know from the studies of LaRu2Si2 f|. The Hund's 
rules ground multiplet is 3 H4 with J = 4, S — 1, L — 5 and the Lande factor 
gi = 4/5. The general Hamiltonian contains the single-ion-like and the intersite 
terms g [7J: 

H = E E KO™(J, J,) + nglnl {-J (J) + i (J) 2 ) + g LftB J ■ B ext (1) 

The first term is the CEF Hamiltonian written for the lowest multiplet given 
by Hund's rules 3 H4. The second term represents the exchange interactions be- 
tween the Pr ions written in the mean-field approximation with the molecular- 
field coefficient n. The third one describes the Zeeman effect. 

3 Results and discussion 

The energy level scheme of the Pr ion in PrRu2Si2 for the tetragonal symmetry 
CEF Hamiltonian contains 5 singlets and 2 doublets ||. A full set of CEF 
parameters relevant to the tetragonal symmetry is given by: B\ = -22 K, B® = 
+ 0.22 K, Bj = + 0.20 K, B% = -12 mK and B\ = -45 mK. This set of parameters 
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has been derived by our self-consistent analysis of experimental data. This set of 
CEF parameters describes well 

i) the energy separations revealed by the inelastic neutron scattering, 

ii) anisotropic temperature dependence of the susceptibility in the paramag- 
netic region, 

iii) the direction of the magnetic moment (along the tetragonal c axis) in the 
ordered state, 

iv) the spontaneous moment of the Pr ion in 4.5 K of about 2.7 \Ib, 

v) the giant magnetocrystalline anisotropy of the magnetization curves at 4.5 

K, 

vi) temperature dependence of the specific heat with the A— type peak at T c . 
The consistent description of so many physical properties provides the strong 

argument for our theoretical approach. 

The ground state is a singlet in the form: 

= 0.703 |+4) + 0.109 |0) + 0.703 |-4) . (2) 

The resulting energy- level scheme is shown in fig.l. It provides excitations 
observable in inelastic-neutron-scattering (INS) experiments at energies of 30 K 
(rg> -> T t2 ) and 330 K (T t2 -> and T t3 — > r^-* , r t4 ). These excitations have 
been observed, indeed, in INS experiments by Mulders et al. |], |j] (2.25 and 
29 meV). In this experiment the 330 K excitation decreases with the increasing 
temperature. This fact we take as the further confirmation of our scheme. With 
the increasing temperature the contribution T t3 — > r t4 , having the lower energy, 
increases together with the Boltzmann population of the T f3 state. This obser- 
vation provides indirect argument for the existence of the state rj^ at 58 K - 
the neutron excitations to it from the lower states are prohibited. As will be 
shown later the existence of a state at about 60 K is in very good agreement with 
the overall temperature dependence of the specific heat. The entropy analysis 
performed in Ref. 1 shows that there is a place for this extra state owing to 
the increasing difference between the experimental and calculated, for two states 
only, entropy as one can see in Fig. 2 of Ref.l. 

Closely lying first excited state together with the ground state creates the 
interesting system. Their field-induced J z components have opposite signs what 
lead to their opposite interactions with the magnetic field. As a consequence 
this 2-levels system behaves as the quasi doublet. Moreover, there is the large 
matrix element between these states. Thanks it the magnetism of this system 
can be relatively easily induced by spin-dependent interactions despite of general 
non-magnetic character of singlets. Such the charge-formed ground state allows 
the appearance of the large magnetic moment, of 2.7 in agreement with 
experimental observation. Also the direction of the moment is related with the 
shape of the eigenfunctions of this 2-level system. 
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3.1 Magnetic susceptibility 



In the quantum paramagnetic theory the magnetic moment is a property of the 
electronic state of the paramagnetic ion. The value of the magnetic moment of 
the ion reflects dependence of the energy of the ion with respect to the magnetic 
field B according to the definition ||: 

E(T) is the total energy of the system over available electronic states resulting 
from the Hamiltonian 1 and shown in Fig.l. The population of states shown 
in Fig. 1 is given by the Boltzmann statistics. The magnetic susceptibility in 
the paramagnetic region was calculated from the direct diagonalization of the 
Hamiltonian, Eq. 1, without the second term. 

The inverse of calculated temperature dependence of the magnetic suscepti- 
bility in the paramagnetic region is shown in Fig.2. It exhibits very anisotropic 
behavior with the easy c axis. The behavior x _1 (T) in the magnetic field parallel 
to the c axis almost follow the Curie- Weiss law. The derived effective param- 
agnetic moment at room temperature amounts to 4.18 /ig. It is larger than the 
free-ion value of 3.58 \xb and this fact is caused by strong CEF interactions. This 
paramagnetic moment decreases with the increase of temperature approaching 
the theoretical value at temperatures comparable with the overall size of the 
CEF interactions. In the perpendicular field the calculated susceptibility is very 
small. The visible discrepancy between calculated and experimental curves we 
attribute to its small value that becomes comparable to the conduction-electron 
contribution. 



3.2 Magnetic state and magnetocrystalline anisotropy 

At T c of 14 K magnetic order wins the temperature disordering. It is accounted 
for in our calculations with the effective exchange interaction parameter n of 2.35 
T/n B . The magnetic ordering produces the abrupt change in the energy level 
scheme that manifests in the A peak at T c . The calculated moment of the Pr 
ion is parallel to the c axis in agreement with experiment. It amounts to 2.7 
Hb at K fully reproducing the experimental datum [[!]]. For completeness we 
add that it turns out from our calculations that at K the Pr 3+ ion in PrRu2Si2 
experiences the internal molecular field of 6.4 T that originates from the intersite 
RKKY interactions. 

The magnetic moment of the Pr 3+ ion is strongly tied to the tetragonal c axis 
as we see from the magnetization curves. In Fig. 3 the full magnetization curves, 
at 4.5 K, calculated within our CEF approach are presented. They are highly 
anisotropic: the field of 5.5 T applied along the c-axis induces (J z ) of 3.6 what 
is about 100 times more then the value for (J x ). This anisotropy preserves also 
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in measurements up to 35 T, where magnetization along the c and a axes is 3.08 
Hb and 0.39 fig , respectively [|TJ. This large anisotropy is in agreement with 
single-crystal magnetic measurements of Shigeoka et al.||. The anisotropy field 
derived from our calculated curves amounts to 400 T. This value corresponds to 
the magnetocrystalline anisotropy energy of 59 J/cm 3 , that is extremaly large. 

For calculations of the magnetization curves in the ferromagnetic state along 
the c-direction, shown in Fig. 3, the primary field-induced moment curve m(Bf u u), 
shown in the insert, has been used. The extraction of the influence of the ex- 
ternal field B ext , the internal field Bi nt = n ■ m has to be subtracted what leads 
to the metastable region at zero field. Such the metastable situation means the 
spontaneous formation of the ferromagnetic state. 

The derived CEF parameters account also for the ferromagnetic ordering with 
its direction, the ordering temperature and the size of the Pr-ion moment and the 
giant magnetic anisotropy. These calculations reveal that the giant anisotropy, 
even above 400 T, can be realized by the crystal-field interactions. 



3.3 Specific heat 

The magnetic component of rare-earth specific heat is calculated by making use 
of the general formula @]: 

cm = <4) 

F(T) is the free energy of the rare-earth subsystem over the available energy 
states resulting from the consideration of the Hamiltonian (eq. 1). The calculated 
contribution of the Pr 3+ ions to the specific heat of PrRu2Si2 is shown in fig.4. 
The overall behavior c±f(T) is in very good agreement with experimentally derived 
magnetic contribution to the specific heat [|I|]. Our calculations reproduce well: 

i) a A-type peak at 14 K associated with the occurrence of the ferromagnetic 
order. There is slight discrepancy in description of the shape of the peak; the 
calculated one is wider then the experimental one. It may indicate more dramatic 
character of the phase transition at T c , 

ii) a tail of the Schottky-like peak above T c . 

The reproduction of the A peak is important outcome of the present paper. 
The previous calculations of Ref. 1 could not get the A peak because they worked 
out with thermally independent electronic states. In fact, the model of Ref. 1 could 
not produce the ferromagnetic order. 

This good agreement indicates on the substantial confidence to the presently 
derived fine electronic structure. Our scheme gives much better description of 
the overall C4/(T) dependence then previous calculations presented in Ref. 1. 
They worked with the temperature independent structure of states as they did 
not take into account the intersite spin-spin interactions and they could not get 
the A peak at T c . 
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4 Conclusions 



The energy level scheme of the Pr 3+ ion in PrRu2Si2 has been constructed on basis 
of all known experimental data. Our scheme is in agreement with the inelastic- 
neutron-scattering data. Our calculations reproduce well the zero-temperature 
moment and its temperature dependence, temperature dependence of the mag- 
netic susceptibility and the specific heat with the sharp peak at T c , single- 
crystalline magnetization curves with the enormous anisotropy field of 400 T. 
The obtained consistent description indicates on the substantial confidence to 
the present evaluation of the crystal field interactions in PrRu 2 Si 2 and should be 
the good starting point for the CEF analysis of other members of RERu 2 Si 2 class 
of compounds. 
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Figure Captions 

Fig. 1. The energy level scheme of the Pr 3+ ion in PrRu 2 Si 2 with the expec- 
tation values of J z and eigenvectors in the paramagnetic state. The tetragonal 
CEF interactions split the 9-fold degenerate ground multiplet 3 H 4 into 5 singlets 
and 2 doublets. Three arrows indicate the allowed INS excitations. They are 
assigned to those experimentally detected in Refs 1 and 4. 
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Fig. 2. Temperature dependence of the magnetic susceptibility of PrRu 2 Si 2 
along the tetragonal c- and a-axis shown in the x^ 1 vs T plot in the paramagnetic 
state. The solid lines show our calculations, points are experimental data after 
Ref. 3. 

Fig. 3. Magnetization curves, the moment vs external field B ext , at 4.5 
K (lines) on single-crystalline PrRu 2 Si 2 calculated along main crystallographic 
directions of the tetragonal unit cell in the ferromagnetic state. Points denote 
experimental data from Ref. 1. In the insert the primary field- induced magnetic 
moment is shown. In the ferromagnetic state, for the extraction of the influence 
of the external field, the internal field B int = n ■ m is subtracted. 

Fig. 4. Temperature dependence of the 4/ specific heat contribution in 
PrRu 2 Si 2 . The dashed-point line shows the result of our calculations. The solid 
line represents calculations of Ref. 1. Points are experimental data from Ref. 1. 
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